Exchange-energy functionals for finite two-dimensional systems 
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Implicit and explicit density functionals for the exchange energy in finite two-dimensional systems 
are developed following the approach of Becke and Roussel [Phys. Rev. A 39, 3761 (1989)]. Excellent 
agreement for the exchange-hole potentials and exchange energies is found when compared with 
the exact-exchange reference data for the two-dimensional uniform electron gas and few-electron 
quantum dots, respectively. Thereby, this work significantly improves the availability of approximate 
density functionals for dealing with electrons in quasi-two-dimensional structures, which have various 
applications in semiconductor nanotechnology. 



I. INTRODUCTION 

Since the advent of density-functional theory (DFT)i 
much effort went in the development of approximate 
functionals for the exchange and correlation energies. 
Most of this work focused on three-dimensional (3D) 
systems, where considerable advances beyond the com- 
monly used local (spin) density approximation (L(S)DA) 
were achieved by generalized gradient approximations, 
orbital functionals, and hybrid functionalsi^ For two- 
dimensional (2D) and low-dimensional systems in gen- 
eral, such efforts have been relatively scarce, partly due 
to the lack of direct applications before the rapid devel- 
opments in semiconductor technology in the early 1980s. 
Presently, the variety of relevant low-dimensional sys- 
tems is very large, including, e.g., modulated semicon- 
ductor layers and surfaces, quantum-Hall systems, spin- 
tronic devices, and quantum dots (QDs). The research 
of electronic, optical, and magnetic properties of low- 
dimensional structures form already a significant contri- 
bution to condensed-matter and materials physics. 

Semiconductor QDs are finite quasi-2D electron sys- 
tems confined typically in GaAs/AlGaAs heterostruc- 
tures.'^ Due to their controllability in size, shape, and 
number of confined electrons, they have an increasing 
number of potential applications in quantum-information 
technology. Within the DFT approach, QDs are most 
commonly treated using the 2D-LSDA exchange func- 
tional derived in 1977 by Rajagopal and Kimball,^ 
which is then combined with the 2D-LSDA correlation 
parametrized first by Tanatar and CeperleyS, and later, 
for the complete range of coUinear spin polarization, 
by Attaccalite and co-workers.^^ Despite the relatively 
good performance of LSDA with respect to, e.g., quan- 
tum Monte Carlo calculations,^ there is lack of accurate 
2D density functionals. Moreover, previous studies have 
shown that functionals developed for 3D systems per- 
form poorly when applied to 2D systems. ^'^'^"^ The exact- 
exchange functional employed within the optimized ef- 
fective potential method, which automatically conforms 
to various dimensionalities, seems an appealing alterna- 



tive, and it has recently been applied to QDs.^^ In that 
method, however, the development of approximations for 
the correlation energies compatible with exact-exchange 
energies remains a complicated problem. 

In this work we present an ab initio derivation of new 
exchange functionals for finite 2D systems in the frame- 
work proposed by Becke and Roussel. We calculate 
the cylindrical average of the exchange hole for a single- 
electron wave function of a 2D harmonic oscillator, and 
use it as a basis of an averaged exchange hole for a generic 
A^-electron system. By utilizing the short-range behavior 
of the exchange hole we are able to derive one implicit 
and one explicit density functional for the exchange en- 
ergy. Both functionals lead to accurate exchange-hole 
potentials and exchange energies when compared with 
the reference data of the 2D electron gas (2DEG) and 
few-electron QDs. The new exchange- energy functionals 
constitute a natural basis for developing corresponding 
correlation-energy functionals. 



II. FORMALISM 

Within the Kohn-Sham (KS) method of spin-DFT,il 
the ground state energy and spin densities p| (r) and 
Pl (r) of a system oi N — N-^ + N i interacting electrons 
are determined. The total energy, which is minimized to 
obtain the ground-state energy, is written as a functional 
of the spin densities (in Hartree atomic units) 

Ev[p],Pi] = Ts[p^,p{\+ E]i[p\+ E,,c[p],Pi\ 

+ I dvv,{v)p„[r), (1) 

where T's[p|, Pi\ is the KS kinetic energy functional, ^^(r) 
is the external (local) scalar potential acting upon the 
interacting system, i?H[p] is the classical electrostatic or 
Hartree energy of the total charge density p(r) — p|(r) -|- 
Pl(r), and Ea:c[p-\,Pi\ is the exchange-correlation energy 
functional. The latter can be further decomposed into 



2 



the exchange and correlation parts as 

Exc[p\,Pi\ = E^[p^^,pi]+ Ec[p],Pi]. (2) 
The exchange- energy functional can be expressed as 



2 ^ 



dvi / dv2 



ri -r2 



-/ij(ri,r2), 



. (3) 

where, within the restriction that the noninteracting 
ground state is nondegenerate and hence takes the form 
of a single Slater determinant, the exchange-hole (or 
Fermi-hole) function ftj is given by 
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(4) 



The sum in the numerator is the one-body spin-density 
matrix of the Slater determinant constructed from the 
KS orbitals, V'fc.o-- The exchange- hole function as defined 
here is always positive. Moreover, integrating this func- 
tion over T2 yields 



j dT2K{Ti,T2) = 1. 



(5) 



This exact property reflects the fact that around an elec- 
tron with spin a at ri , other electrons of the same spin 
are less likely to be found as a consequence of the Pauli 
principle. From Eq. ([3]) it is clear that to evaluate the 
exchange energy in 2D, we just need to know the cylin- 
drical average w.r.t. s = r2 — ri of the exchange-hole 
around ri. Expressing the exchange-hole by its Taylor 
expansion 



K{vi,V2 = n +s) = exp(s • V')/i"(ri,r')|r'=ri, 
the cylindrical average is defined as 
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.exp(s.V')/^J(ri,r')|. 



(6) 
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The short-range behavior with respect to s is then ob- 
tained as 



= p,{v,) + s^C:{v,) + ...^ 



ri + 



(8) 



where is the so-called local curvature of the exchange 
hole around the given reference point ri. This function 
can be expressed a o^^'^'^i^^ 



where 
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1 (Vp.) 

2 Pa 



(9) 
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is (twice) the spin-dependent kinetic-energy density, and 

jp.- = ^ E (^^'».-) - (^^l-) (11) 

* k=\ 

is the spin-dependent paramagnetic current density. 
Both To- and jp^o- depend explicitly on the KS orbitals. 
Thus the expression in Eq. ([9]) has an implicit dependence 
on the spin-densities pa- Now, once we can provide an 
approximation for /ij(ri,s) satisfying the normalization 
condition of Eq. ([S]) , we can compute the exchange-hole 
potential 

/"OO 

C/,^(ri) = -2^ / dsKiv^.s), (12) 
Jo 

and finally the exchange energy 

EAP^,Pi\ ^ ^E/ rfrip.(ri)[/;(ri). (13) 

III. EXCHANGE-ENERGY FUNCTIONAL 
A. Model 

As the basis of our exchange model we consider the 
ground state single-electron wave function of a 2D har- 
monic oscillator 



V'a-(r) = —j= exp 



(14) 



which is normalized for each a 7^ 0. We point out that 
the harmonic (parabolic) approximation for the confin- 
ing potential is the most common choice when modeling 
QDs fabricated in semiconductor heterostructuresi^ The 
exact exchange-hole function Q for the single-particle 
case becomes /i^(ri,r2) = ■\\)*„(y:2)i}cr(y2) = Pa-(r2)- Set- 
ting ri = r and r2 = r + s, we calculate the cylindrical 
average as 



1 f^- 

TT- \ d(l)sPa{r + s) 
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X / d0s exp [— 2a^rscos((/))] 



a 



■ exp 



-a^ (r^ + s2)] Ioi2a^rs), (15) 



where Io{x) is the zeroth order modified Bessel function 
of the first kind (note that /o(0) = 1). Performing the 
integral in Eq. ([T^ yields 

U^ir) = -|a|V^exp [-aVV2] /o(aVV2) (16) 

for the exchange-hole potential. Since the modi- 
fied Bessel function has the limit property Io{x) — s- 
exp(a;)/-\/27rx, when x — > +00 with a; G M, we imme- 
diately find lim^^oo U^{r) — — l/r. 



B. Implicit density functional 



At this point, we adopt the strategy of Becke and Rous- 
seli^ and elevate expression (fT5)) as a general model for 
the averaged exchange hole of a generic A^-electron 2D 
system. In order to locally reproduce the short-range 
behavior of the exchange hole, we replace o? and by 
functions of r, respectively, i.e., o? a(r) and H^)- 
Now we can rewrite Eq. (fT5|) as 



Kia,b;s) = ^exp [-^a{b + s^)] h){2aVhs). (17) 

This model satisfies, through its original definitions, 
h'^{a,b;s) > and the unit normalization condition of 
Eq. ([5]). From the second-order term in the Taylor ex- 
pansion in Eq. ([8]) we obtain 



(y - l)exp(?/) 



where y := ab. The zeroth-order term gives 
a = Tipa exp(y). 



and hence we get 



exp(-?/). 



(18) 



(19) 



(20) 



As the result, Eqs. (|T8)) - ([20|) determine, together with 
Eqs. P?|l . (fT^ and an implicit density functional. 

This completes the derivation of our first approximation 
to the exchange-energy functional. 
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FIG. 1: Exchange hole (a) and the exchange energy per par- 
ticle (b) for two-dimensional electron gas (^ = 0) in the exact 
expression (dashed line) and in the present model (solid line). 



This expression can be compared with the exact 2DEG 
result widely applied to finite systems in terms of the 
LSDA, 



,2DEG 
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37rr., 



(l+0'^2 + (l-0^/2l . (23) 



Interestingly, the only difference between these two ex- 
pressions for Cx is in the prefactor, which is about ~ 4.4% 
smaller in the latter. Concluding, Eq. (p2)) can be used as 
an explicit density functional in the LSDA fashion when 
applied to inhomogeneous systems. 



C. Explicit density functional 



We now show that one can also derive an explicit 
density functional. For that purpose, we consider the 
2DEG where the derivatives in Eq. ^ are zero, i.e., 
V^Pct = and Vpcr — 0, and we take the known ex- 
pression for the kinetic- energy density of the 2DEG,iS, 
Tcr = 27rp^ -I- jp^/po-- Note that the current-dependent 



term in Eq. ([9]) cancels the current-contribution to the 
kinetic energy, leading to the expected result that the 
exchange energy of a uniform electron gas does not de- 
pend on the current density^^ These simplifications lead 
to y = 0, 6 = 0, and a = irpa-, so that the averaged 
exchange-hole function becomes 



(21) 



and the exchange-hole potential is now given by U^[pa\ = 

1/2 

—'KpJ , which is an explicit functional of the spin- 
density. Defining the 2D density parameter — l/y^Trp, 
where p = p^ + pi is the total density, and the polar- 
ization ^ = [p] — Pi)/p, the total exchange energy per 
particle becomes 
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(l + C)'/' + (l-0'/' 



(22) 



IV. NUMERICAL RESULTS 

A. Two-dimensional electron gas 

Besides the exchange energies, it is interesting to com- 
pare the averaged exchange hole given in Eq. (j2ip with 
the exact exchange hole of the 2DEG. Following the 
derivation of Gori-Giorgi and co-workers^i for the pair- 
distribution functions of the 2DEG, we find 



l2DEG 



(s) = 4p, 



(24) 



where klp^ 



(47rpo.)^/'^ is the Fermi momentum (for spin 
(t) in 2D, and Ji is the ordinary Bessel function of the first 
order. In Fig.[T](a) we compare the exchange holes with a 
fixed spin-density po- = 1 between our model (solid line) 
and the exact 2DEG result (dashed line) and find a good 
qualitative agreement. Figure [T](b) demonstrates the dif- 
ferences in the exchange energies (per particle) given in 
Eqs. (|22p and (|23p as a function of the total density (for 
C = 0). 
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FIG. 2: Exchange-hole potential of a two-electron quantum 
dot calculated exactly (dashed line) and using the implicit 
density functional of Sec. IIIIBI (solid line) , standard LSDA 
(dotted line), and the explicit density functional of Sec. IIII Cl 
(dash-dotted line) . The analytic density and orbitals are used 
as input for the model and LSDA. 



B. Few-electron quantum dots 

Next, we consider the smallest nontrivial QD consist- 
ing of two electrons. The Hamiltonian is given by 



1 



r2 - ri| 



(25) 



where we set the strength of the harmonic confinement 
to (jjQ ~ 1. In this case the ground-state (singlet) wave 
function is known analyticallyji-^ and the total density 
can be expressed as^ 



p{r) 



7r(V27r + 3) 
/( 



{e-'-'(l + rV2) + 



1 



/o(rV2) + r2/o(rV2) + r^h{r''/2)\ }. (26) 
The exact (spin) exchange-hole potential is simply 



r - r' 



(27) 



In Fig. [2] we compare the exact exchange-hole potential 
(dashed line) to the result of our implicit density func- 
tional (solid line) and find an excellent agreement. We 
note that in the range 0.75 ^ r < 0.85 the parameter 
y is not solvable from Eq. (HH]). Therefore, we set y to 
zero in this region. This corresponds to the 2DEG result 
discussed above, i.e., we perform a well- valid LSD A- type 
approximation in this small regime. The dash-dotted line 
shows the result of our explicit density functional which 
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FIG. 3: Similar to Fig. [2] but for 20 electrons. The exact- 
exchange result was calculated using the KLI approxima- 
tion. Self-consistent (correlation included) density and or- 
bitals from the standard LSDA are used as input for our 
functionals and LSDA exchange. 



TABLE I: Exchange energies for different number of electrons 
calculated using the exact exchange, the implicit density func- 
tional of Sec. IIII 51 standard LSDA, and the the explicit den- 
sity functional of Sec. IIII CI For A'' > 2 the EXX result was 
calculated within the KLI approximation, and as input for our 
functionals and LSDA exchange we used the self-consistent 
density and orbitals from the standard LSDA. 



N 


EXX 


Implicit functional 


LSDA 


Explicit functional 


2 


-1.0839 


-1.0836 


-0.983 


-1.026 


6 


-2.229 


-2.284 


-2.130 


-2.223 


12 


-4.890 


-5.059 


-4.763 


-4.972 


20 


-8.781 


-9.124 


-8.632 


-9.012 



is identical to the implicit functional in the 2DEG limit. 
This density functional has the wrong asymptotic behav- 
ior, but it is considerably closer to the exact result than 
that of the standard LSDA of Eq. ^ (dotted hue). 

In Fig. [3] we show the exchange-hole potentials for a 
20-electron QD with = 0.42168, corresponding to a 
typical confinement of 5 meV in the effective-mass ap- 
proximation when modeling QDs in GaAsi^ The exact- 
exchange (EXX) result (dashed line) is calculated in 
the Krieger-Li-Iafrate (KLI) approach'^*' within spin-DFT 
implemented in the octopus real-space code.^^ Again, we 
find a very good agreement between the EXX and the re- 
sults of our implicit density functional (solid line). Also, 
the explicit functional (dash-dotted line) reproduces the 
EXX curve at small r by a good accuracy, whereas the 
LSDA curve (dotted line) lies considerably above the 
EXX result in this regime. However, as seen in Table |I] 
the LSDA exchange energy is rather accurate. This is 
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due to the fact that the difference in the EXX and LSDA 
curves at small r is compensated at 5 < r < 6, where 
the LSDA curve lies below the EXX one. In smaller 
QDs instead, the table shows that our exchange-energy 
functionals are superior to the LSDA also in terms of the 
exchange energies, i.e., not only in terms of the exchange- 
hole potentials. 

V. CONCLUSIONS 

To conclude, we have provided new exchange function- 
als for finite two-dimensional systems that, until now, 
have lacked accurate functionals beyond the local den- 
sity approximations. We have presented ah initio deriva- 
tion of both implicit and explicit exchange functionals 
that significantly improve the exchange-hole potentials 



and exchange energies of few-electron two-dimensional 
quantum dots. Our results also suggest that accurate 
approximations for the correlation energy functional can 
be developed analogously. 
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